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ABSTRACT 



We numerically investigate abundance properties of the Galactic globular clusters 
(GCs) by adopting a new "external pollution" scenario. In this framework, GCs are 
assumed to originate in forming low-mass dwarfs embedded in dark matter subhalos 
at very high redshifts [z) and thus be chemically influenced by field AGB stars of the 
dwarfs during early GC formation processes. GCs within a dwarf galaxy therefore can 
' be formed from the mixture of (i) gas ejected from the field AGB stars formed earlier 

, in the dwarf and (ii) the interstellar gas infalling to the central region of the dwarf. 

^SJ ■ In this external pollution scenario, the ratio of the total mass of infalling gas to that 

I of AGB ejecta during GC formation in a dwarf (s) and the time scale of gas infall 

. (iTi) are the most important key parameters that can determine abundance properties 

' of GCs. We mainly investigate the abundance inhomogeneity among light elements 

fH \ (e-g-) C, N, O, Na, and Al) of stars in GCs by using the latest stellar yield models of 

metal-poor AGB stars with and without third dredge-up. Our principal results for the 
I ' models with no third dredge-up, which are more consistent with observations, are as 

? , follows. 

(1) Both [N/Fe] and [C/Fe] can be diverse among stars within a GC owing to 
chemical pollution from field AGB stars. [N/Fe] distributions in some GCs can clearly 
show bimodality whereas [C/Fe] is monomodal in most models. [N/Fe] distributions 
depend on s such that models with smaller s (i.e., larger mass fraction of AGB ejecta 
used for GC formation) show the [N/Fe] bimodality more clearly. 

(2) N-rich, C-poor stars in GCs also have higher He abundances owing to pollution 
' from massive AGB stars with He-rich ejecta. The number fraction of He-rich stars (Y 

> 0.30) is higher for the models with smaller s and shorter crj for 3 < s < 24 and 
10^ < cTi < 10^ yr. He abundances of stars correlate with [N/Fe] and [Al/Fe] and 
anticorrelate with [C/Fe], [0/Fe], and [Na/Fe] within GCs in our models. 

(3) Although our model can much better explain the observed C-N and Mg-Al 
ant icor relations than previous theoretical models, it is in strong disagreement with 
the observed 0-Na anticorrelation. 

(4) This model naturally provides an explanation for the large fraction of CN- 
strong stars without recourse to an implausible IMF. 

Based on these results for the above external pollution scenario, we discuss the 
long-standing problem of the CN-bimodality prevalent in the Galactic GCs, the pos- 
sible helium abundance inhomogeneity in these systems, and their horizontal branch 
morphologies. 

Key words: globular clusters: general - globular clusters: individual {uj Centauri)- 
globular clusters:individual (NGC 6752)-- galaxies: star clusters - galaxies: evolution - 
galaxies: stellar content 



1 INTRODUCTION 
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Since observational evidence of star-to-star abundance inho- 
mogeneity among the light elements of stars in the Galactic 
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Table 1. Stellar yields of AGB stars 
model H He C N O Na 24]y[g 25 ^jg 26y^^ ^1 Pg 

YAl .6639E+0 .3345E+0 .5316E-4 .1143E-2 .1592E-4 .4365E-6 .2806E-6 .2786E-4 .3429E-5 .1855E-4 .5481E-4 
YA2 .6414E+0 .3491E+0 .1206E-2 .7358E-2 .2475E-3 .1438E-4 .9788E-6 .8031E-4 .1121E-3 .1067E-4 .5375E-4 
Notes: The mass fraction of AGB ejecta is given for each element. These yields are for AGB stars with masses of 6.5M0 and [Fe/H] = 
— 1.4. Stellar yields for the AGB stars with differences masses and metallicities are from Campbell et al. (2006). 



Table 2. Initial abundances of infalling gas 



model 


Y 


[C/Fe] 


[N/Fe] 


[O/Fc] 


[Na/Fe] 


[Mg/Fe] 


[Al/Fo] 


YIl 


0.23 


0.0 


-0.25 


0.3 


-0.3 


0.5 


0.1 


YI2 


0.23 


1.25 


-1.97 


-1.63 


-0.3 


0.38 


-1.74 


YI3 


0.23 


0.0 


-0.25 


0.3 


-0.3 


0.5 


0.3 


YI4 


0.23 


-0.3 


-0.25 


0.3 


-0.3 


0.5 


0.1 



globular clusters (GCs) was discovered (e.g., Cohen 1978; 
Peterson 1980; Norris et al. 1981; Loop et al. 1986), the origin 
of the inhornogcnoity has been extensively discussed both 
theoretically and observationally (e.g., Sneden et al. 1992; 
Norris & Da Costa 1995; Smith et al. 2005; see Gratton et al. 
2004 for a recent review) . So far the following two scenarios 
(or working hypotheses) have been proposed for the origin 
of the abundance inhomogeneity: (1) the primordial hypoth- 
esis and (2) the mixing hypothesis (e.g., Sweigart & Mengel 
1979; Cottrell & Da Costa 1981; Freeman & Norris 1981; 
Smith 1987; Suntzeff 1993; Kraft 1994; Denissenkov & Weiss 
1996). The first hypothesis is that the observed inhomogene- 
ity is due to the second generation of stars that formed from 
gas ejected from the first generation of evolved stars (e.g., 
AGB stars) ("primordial pollution" or "self-pollution" sce- 
narios, e.g., Cottrell & Da Costa 1981). 

Internal processes of stars, such as dredge-up of the CN- 
processing of envelope material in inner hydrogen-burning 
regions, are key for understanding the observed chemical in- 
homogeneity of GCs in the second mixing hypothesis (e.g., 
"deep mixing"; Smith 1987; Kraft 1994; Thoul et al. 2002). 
Although deep mixing is highly unlikely to occur in less 
evolved staxs on the main sequence and subgiant-branch, 
star-to-star abundance inhomogeneity was found in these 
stars of some Galactic GCs (e.g.. Cannon et al. 1998 for 47 
Tuc). The self-pollution scenario is accordingly suggested to 
be more promising than the mixing (evolutionary) scenario 
(Da Costa et al. 2004; Gratton 2004 for a recent review). 

Stellar ejecta from AGB stars and massive staxs have 
been considered to play a key role in early chemical evo- 
lution of GCs for the self-pollution scenario (e.g., Fenner et 
al. 2004; Charbonnel & Prantzos 2006; Prantzos & Charbon- 
nel 2006; Karakas et al. 2006; Bekki & Chiba 2007). Previ- 
ous studies suggested that GCs are unlikely to retain AGB 
ejecta owing to ram pressure stripping by warm/cold gas of 
the Galactic halo and disk (e.g., Frank & Gisler 1976; Smith 
1996; Gnedin et al. 2002). These imply that GCs are even 
more unlikely to retain AGB ejecta in protogalactic environ- 
ments, where denser halo/disk gas can strip the AGB ejecta 
efficiently from GCs during hydrodynamical interaction be- 
tween the gas and the AGB ejecta. These previous studies 
thus appear to suggest that self-pollution is not likely within 
GCs evolving in isolation, though numerical attempts have 
not yet been made to investigate the details of self-pollution 
processes within GCs. 

Furthermore, the self-pollution scenario is suggested to 
have the following two problems in explaining physical prop- 



erties of GCs. The first is that the observed large number 
fraction of CN-strong populations (^ 0.5 for NGC6752), 
which correspond to the second generation of stars formed 
from AGB ejecta, can not be explained by the scenario with- 
out assuming very unusual and unrealistic IMFs (e.g.. Smith 
& Norris 1982; D'Antona and Caloi 2004): the total mass of 
AGB ejecta from the first generation of stars in a GC is 
only 1 — 10% of the total mass of the GC with a canoni- 
cal IMF (Bekki & Norris 2006). The second is that the ob- 
served O— Na and Mg— Al anticorrelations and C-depleted 
stars in CN-strong populations can not be reproduced quan- 
titatively by chemical evolution models based on the self- 
pollution scenario (Fenner et al. 2004). This problem, how- 
ever, could be due to the adopted AGB models, which have 
some uncertainties in predicting chemical yields (e.g., Denis- 
senkov & Herwig 2003; Karakas & Lattanzio 2003; Ventura 
& D'Antona 2005a,b; Campbell et al. 2006). 

The most serious problem among these is the first one 
- the observed large fraction of stars that (in this scenario) 
formed from AGB ejecta. The very large mass fraction of 
"polluted" second generation of stars in the self-pollution 
scenario requires that AGB ejecta used for the formation of 
this generation of stars originate from stellar systems that 
had total masses much larger than those of the present GCs 
and have already disappeared for some reason. There could 
be two possible scenarios that can satisfy the above require- 
ment. In the first place, we propose: (i) the GCs were ini- 
tially more massive than the present-day GCs by a factor of 
10— 100, (ii) AGB ejecta of the first generation of stars were 
consumed by the formation of the second generation with 
very high star formation efficiency, and (iii) most (90 — 99%) 
of the first generation of stars had been preferentially lost by 
some physical mechanisms such as tidal stripping due to the 
Galactic tidal field. This scenario can be viable if enough of 
the AGB ejecta can be converted into the second generation 
in less than 10* yr after the completion of the first genera^ 
tion Type II supernovae (SNell) - otherwise, age differences 
between the first and the second generation would cause a 
large age spread inconsistent with the observed tightness of 
the color-magnitude diagrams of GCs. 

In the second scenario, proposed by Bekki (2006): (i) 
proto-globular clouds are located in the central regions of 
low-mass proto-galaxies embedded in dark matter subha- 
los, (ii) AGB ejecta of the host galaxy's field stars, which 
formed earlier and surround the clouds, goes into star for- 
mation within the clouds on a time scale of ~ 10^ yr (i.e., 
until the intracluster gas is expelled by SNell) , and (iii) the 
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Table 3. Model parameters 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Model no. 


Gas yield 


AGB yield s (Min/Magb) 


(7i (in units of lO" yr) Mg(0) (in units of Mf) 


Comments 


Ml 


YIl 


YAl 


3.0 


0.1 


0.01 


The fiducial model 


M2 


YI2 


YAl 


3.0 


0.1 


0.01 




M3 


YI3 


YAl 


3.0 


0.1 


0.01 




M4 


YIl 


YA2 


3.0 


0.1 


0.01 


with third dredge-up 


M5 


YIl 


YAl 


6.0 


0.1 


0.01 




M6 


YIl 


YAl 


12.0 


0.1 


0.01 




M7 


Yll 


YAl 


24.0 


0.1 


0.01 




M8 


Yll 


YAl 


3.0 


0.01 


0.01 




M9 


Yll 


YAl 


3.0 


0.05 


0.01 




MIO 


Yll 


YAl 


3.0 


1.0 


0.01 




Mil 


YIl 


YAl 


3.0 


0.1 


0.001 


A model for NGC 6752 


M12 


YIl 


YAl 


3.0 


0.1 


0.005 




M13 


YIl 


YAl 


3.0 


0.1 


0.029 




M14 


YIl 


YAl 


0.17 


0.01 


0.001 


A model for oj Cen 


M15 


YI4 


YAl 


3.0 


0.1 


0.01 


smaller initial [C/Fe] 


Notes: Cols. 


(2-3) yield types for infalling ga 


s and AGB ejecta, (shown in 


tables 1 and 2) respectively. Cols. (4-6) model parameters. 



low-mass galaxies were destroyed by the Galactic tidal fields 
to become the Galactic halo components whereas the GCs 
were tidally stripped to form the Galactic halo GCs. In this 
scenario there can be very small age differences (~ lO'^yr) 
between stars with different abundances in a single GC. This 
pollution process would be better called "external pollution" 
rather than "self-pollution", because most of AGB ejecta 
originates from outside the GCs. Although Bekki & Norris 
(2006) have recently discussed the origin of the possible he- 
lium overabundance of stellar populations of co Cen in the 
context of the external pollution scenario, the abundance 
inhomogeneities observed in "normal" GCs have not been 
extensively studied using previous chemical evolution mod- 
els based on the external pollution scenario. 

The purpose of this paper is thus to discuss the origin of 
abundance inhomogeneity observed in the lighter elements 
(e.g., C, N, O, Na, and Al) of GCs by using chemical evolu- 
tion models based on the external pollution scenario. In par- 
ticular, we discuss (1) the [C/Fc]-[N/Fe] anticorrelation, (2) 
CN-bimodality, (3) the O— Na and Mg— AI anticorrelations, 
(4) the helium abundance, and (5) observable predictions of 
the models. We adopt the AGB models recently developed 
by Campbell et al. (2006) in which the third dredge up in 
low-mass AGB stars can, in efi^ect, be switched ofi^ and on 
due to the use or non-use of the Ledoux criterion for convec- 
tive boundaries. We mainly show the results of the models 
without the third dredge-up, because the models can better 
explain the observed abundance inhomogeneity if they are 
included in the present chemical evolution models. We re- 
alize that these models are ad-hoc, but they enable us to 
investigate the proposed formation mechanism. 



2 THE MODEL 

2.1 The external pollution scenario 

We consider GC formation in the central regions of low- 
mass, proto-galaxies embedded in dark matter halos at very 
high z. The GC host galaxies can provide a deep gravita- 
tional potential that allows AGB ejecta to be retained and 
used for star formation in proto-GC clouds. In this scenario, 
both (i) AGB ejecta from field stellar populations (i.e., the 



main components of the galaxies) and (ii) proto-galactic in- 
terstellar gas (which falls into the central regions) can be 
mixed with each other and subsequently used for GC forma- 
tion. Gas ejected from SNell with velocities of ~ 1000 km 
s~^ is assumed to be expelled from dwarfs during galaxy for- 
mation and thus irresponsible for the chemical evolution of 
forming GCs. After GC formation, the host galaxies merge 
with the proto-Galaxy and axe completely destroyed by the 
strong Galactic tidal field. The field stellar components of 
the galaxies are dispersed into the Galactic halo region to 
become old, metal-poor halo stars whereas GCs, which are 
not destroyed by the Galaxy owing to their initial compact- 
ness, become the Galactic halo GCs. 

Star formation and chemical evolution are assumed to 
proceed until type II supcrnovac (SNe II) formed within 
proto-GCs prevent further star formation in this scenario. 
Therefore, star formation is assumed to continue for less 
than ~ lO'^ yr in the models in which stars with different 
masses form in a coeval way. Age differences of stars are typ- 
ically less than 5 Myr for most models in the present study. 
Evolution of chemical abundances depends mainly on (1) 
how much AGB ejecta is mixed into fresh ISM and conse- 
quently converted into stars and (2) the time scale of gas 
infall into central regions of galaxies. Advantages and dis- 
advantages of the external pollution scenario are discussed 
later in §4.1. 

2.2 Star formation and gas mass evolution 

The total gas mass (Mg(t)) of a forming GC changes through 
the gas ejection of AGB stars, gas infall, and star formation: 

^ = A{t)+m-s{t), (1) 

where A{t), I{t), and S(t) arc the gas ejection rate of the 
field AGB stars, gas infall rate, and star formation rate, 
respectively. In order to calculate A{t) at each time step for 
the field stars (which have a total mass of Mi), we assume an 
IMF in number defined as tpirni) = Aomi^', where mi is the 
initial mass of each individual star and the slope q = 2.35 
corresponds to the Salpeter IMF. The normalization factor 
Aq is a function of Mf, m\ (lower mass cut-off), and 
(upper mass cut-off): 
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Figure 1. The upper panel shows the time evolution of the gas 
infall rate (/(t), solid) and the gas ejection rate of AGB stars 
{A{t), dotted) in the fiducial model Ml. The lower panel shows 
the time evolution of the total mass of AGB ejecta accumulated 
in the proto-GC cloud (Magb> solid), the gas accreted in the 
cloud via gas infall (Min, dotted), and the gas in the cloud (Mg, 
dashed) for Ml without star formation. The total mass of the 
proto-GC cloud increases both by accretion of AGB ejecta and 
by gas infall. 



^0 = 



Mf X (2 - q) 



2-9 . 



mi 



2-5 • 



(2) 



where mi and ni^ arc set to be O.IM0 and 120Mq, respec- 
tively. We adopt q = 2.35 for all models in the present study. 
A{t) between t and t + dt \s accordingly described as: 



feii>{m)dm, 



(3) 



where mi and m2 are masses of stars that turn off the main 
sequence (thus becomes a AGB star) at the time t + dt and 
t, respectively (thus mi < m2). fej describes the total gas 
mass ejected from an AGB star with initial mass mi and final 
mass (mp). We derive an analytic form of /ej (= mi — mp) 
from the observational data by Weidemann (2000) by using 
the least-square fitting method, and find: 



fej = O.9I6M1 - 0.444. 



(4) 



In order to calculate the main-sequence turn-off mass 
(mxo = mi and m2) at each time step, we use the following 
formula (Renzini & Buzzoni 1986; Norman & Scoville 1988): 



logmTo(ts) = 0.0558(logts) - 1.338 logts + 7.764, 



(5) 



where mro is in solar units and time in years. We assume 
that the time t = corresponds to the epoch (To) when the 
most massive AGB star (with mi = BM©) starts to eject gas 
in the present models. Therefore ts = t + To in the above 
equation (5). 

The gas infall rate I{t) is described as 



I{t) = lo exp( 



2<Tl2 



(6) 



where Jo is the normalization factor for the total mass of 
infalling gas, Ti represents the epoch of the maximum in- 
fall rate, and ai describes the time-scale of gas infall. The 
adopted Gaussian distribution means that the gas infall rate 
increases, whether slowly or rapidly, - depending on the 
adopted CTi, for < t < Ti - and then decreases for Ti < t. 
The star formation rate S{t) is described as 



S{t) = SoMg 



(7) 



where So (a normalization constant) is determined for a 
given a so that most of the gas can be consumed within 
an order of ^ 10^ yr. Although we investigate a range of a, 
we show the results for the models with a = 1 which can 
better explain observations. 

Observations (e.g., tightness of color- magnitude dia- 
grams of GCs) strongly suggest that the timescale of GC 
formation should be quite short (an order of 10* yr). We 
also consider that star formation is truncated owing to feed- 
back effects of Type-II supornovao (SNc II) formed during 
GC formation: the present models with this truncation can 
keep abundance homogeneity in heavy elements (e.g., Fe) 
and thus can be consistent with physical properties of nor- 
mal GCs in the Galaxy. Therefore we assume that gas can 
be converted into stars in GCs only for < t < Tend, where 
Tend should be similar to or shorter than 10^ yr correspond- 
ing to the typical lifetimes of stars that explode as SNe II. 
Wc present the results of the representative models with 
Tend ~ 10^ yr and Tj/Tcnd = 0.5 in the present study. The 
gas mass (Mg) is given in units of Mf (the total mass of field 
stars). In the adopted IMF, Magb is 0.01 (xMf) att = Tend 
for all models. 



2.3 Abundance evolution 

The time evolution of the abundance of the i-th element is 
described as follows: 
d(ZiMg) 



dt 



^A,«(t)A(t) + Zi4t)i{t) - z,s{t), 



(8) 



where ZA,iit) and Zi^iit) are the i-th abundance of AGB 
ejecta and infalling gas, respectively. Z\^i{t) is time- 
dependent, because of the varying timescales of evolution 
for the different masses of AGB stars (see equation (5)). We 
assume that Zi^i(t) is constant for each i-th element and 
thus refer to it simply as Zi^i in the present study. 

In order to calculate ZA,i{t) , we use the latest chemical 
yield tables produced by Campbell et al. (2006) in which 
H, He, C, N, O, Na, ^''Mg, ^^Mg, ^'^Mg, Al, are listed for 
AGB stars with the masses ranging from 2.5M0 to 6.5Mq 
for a given [Fe/H] (—1.4). These models follow normal AGB 
evolution with hot bottom burning, but we produce two sets 
of yields - one with the normal third dredge-up and one 
without (see below for details). The chemical yield model 
without (with) the third dredge-up are referred to as YAl 
(YA2) in the present study. Table 1 shows the mass fraction 
of the above eleven species in AGB ejecta from Campbell 
et al. (2006). Wc consider that Zi^i should be similar to 
typical values of old halo stars in the Galgixy (e.g., Goswami 
& Prantzos 2000; Venn et al. 2004) and the adopted values 
of Zi^i are summarized in Table 2 for four different models. 
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The details of the models without third dredge-up are given 
in Campbell (2006). Throughout this study, we assume that 
[Fe/H] is the same between infalling gas and AGB ejecta. 
This assumption is reasonable, because we consider that gas 
from SNcII can be expelled from dwarf efficiently during 
galaxy formation. We do not include the effects of massive 
rotating stars on the chemical evolution of proto-GC clouds 
in the present study, because we intend to understand more 
clearly the roles of AGB stars in the chemical evolution: we 
need to derive separately the roles of AGB stars and those 
of massive rotating stars in the chemical evolution. 

2.4 With and without third dredge-up 

We mainly show the results of the models that use the 
AGB yields with no third dredge-up, though we investi- 
gate the models with and without third dredge-up. It should 
be stressed here that the models with no third drcdgo-up 
are purely hypothetical ones that have not yet been proved. 
Previous studies based on self-pollution models with third 
dredge-up in AGB stars failed to explain many of the obser- 
vations of abundance inhomogeneity in GCs (e.g., Fenner et 
al. 2004). If our models can successfully reproduce some ob- 
servations that were not explained by previous ones, we can 
obtain some insights into the "culprits" in the self-pollution 
or external pollution scenarios. The abundance inhomogene- 
ity of GCs in the present models with no third dredge-up 
originates from pure hydrogen burning, which is available in 
HBB AGB stars. 

The present models with third dredge-up can not ex- 
plain at all most of the observed correlations between dif- 
ferent abundances (e.g., the [C/Fo]-[N/Fe] anticorrelation) . 
For example, these models can not produce stars with [C/Fe] 
< and [N/Fe] > 0, because [N/Fe] increases with the in- 
crease of [C/Fe] during the chemical evolution owing to the 
C-rich ejecta of AGB stars. Furthermore, the observed CN- 
bimodality of stars can not be reproduced in these models, 
owing to the lack of C-depleted, N-rich stars in these mod- 
els. These failures would suggest that AGB stars are not 
the polluters that are responsible for the observed abundance 
inhomogeneity in GCs. Wc however consider that it is still 
worthwhile to investigate abundance properties of GCs by 
using diflferent AGB models (i.e., those with no third dredge- 
up) in the present study for the above-mentioned reason. It 
should also be stressed here that the present models with no 
third dredge-up can better explain observations, because we 
adopt the star formation histories of proto-GC clouds in the 
external pollution scenario: successful reproduction of some 
observations in the present study is due to the combination 
of the adopted models both with no third dredge-up and star 
formation histories of the, external pollution scenario. 

Further, at present there are no yields available for the 
so-called "Super- AGB Stars", which are more massive AGB 
stars that experience thermal pulses after core carbon burn- 
ing. The few calculations in existence (e.g., Ritossa et al 
1996; Gil-Pons et al 2005; Siess 2006) show that these stars 
experience very hot HBB and very little dredge-up. Quali- 
tatively this is what is required to solve the GC problem, so 
our models without dredge-up may simulate to some extent 
the yields from these stars, while we await full evolution and 
nucleosynthesis calculations for the appropriate masses and 
compositions. 



We close this section with some comments about the 
reliability of AGB star yields. We believe that the balance 
of evidence indicates that AGB stars with dredge-up are not 
the site of the abundance inhomogeneities seen in globular 
cluster stars. However, there are still so many uncertain- 
ties that it may still be possible for these stars to produce 
the abundance profiles required. For example, the depth of 
dredge-up is still very uncertain, the temperature profile in 
the HBB envelope depends on the convection theory (e.g., 
Ventura & D'Antona 2005a, b, c) as well as associated pa- 
rameters, mass-loss and of course various reaction rates. Al- 
though previous studies (Ventura et al. 2002; Ventura & 
D'Antona 2005c) provided useful tables for stellar yields of 
AGB stars, the lack of tables for some abundances (e.g., the 
lack of a table for Al abundance in Ventura & D'Antona 
2005c) currently does not allow us to investigate the 0-Na 
and Mg-Al anticorrelations in a fully self-consistent manner. 
It is thus our future study to compare the present results 
with those based on chemical evolution models using tables 
from works other than Campbell (2006). More comments on 
AGB models with and without third dredge-up are given in 
the Appendix B. 

2.5 Main points of analysis 

One of the most important parameters that govern the 
chemical evolution of GCs in the external pollution scenario 
is the ratio (s) of the total mass of the infalling gas {Mm) 
to that of AGB ejecta (Magb). We investigate the models 
with s (Min/Magb) ranging from 0.1 to 100 and thereby 
isolated the best values of s that can explain observations. 
The time scale of gas infall (ai) is also an important param- 
eter that can control final abundance distributions of GCs 
in the present study. We investigate models with ai rang- 
ing from 10^ yr to lO'^ yr and show important results of the 
models. 

Although these s and ai are the most important param- 
eters, the initial gas mass Mg(0) can also change abundance 
distributions in a less dramatic maimer for a given s and 
cri. We investigate the models with 0.001 < Mg(0) < 0.029 
and show the results of some representative models. Table 
3 summarizes the parameter values used for the 15 models 
investigated in the present study. We first present the results 
of the fiducial model Ml, which shows interesting and im- 
portant results, in §3.1. We then describe the dependences 
of the results on the two key parameters {s and ai) in §3.2. 

By changing s, ai, and Mg(0), we can construct a model 
that can explain reasonably well the abundance distribution 
for a specific GC. For example, the model Mil shows a bi- 
modal [N/Fe] distribution with the normalized number frac- 
tion of N-rich and C-depleted stars being equal to that of 
the N-normal and C-normal ones and thus is consistent with 
the observed CN distribution in NGC 6752 (e.g.. Smith & 
Norris 1993). Since wc are focusing on the general behaviors 
resulting from parameter dependences of the models, we do 
not discuss extensively the abundance distributions of indi- 
vidual GCs. We however just briefiy discuss the origin of He 
abundance distributions of different stellar populations in uj 
Cen in the Appendix A, because this issue is quite topical 
and is currently being investigated extensively by different 
theoretical models (e.g., Bekki & Norris 2006; Maeder & 
Meynet 2006). 
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Figure 2. Time evolution of the fiducial model Ml (solid) and 
M2 (dotted) on the [C/Fe]-[N/Fe] plane. The initial and final 
values are plotted by open and filled squares, respectively. For 
comparison, observational results by Smith &; Norris (1993) for 
AGB stars in NGC 6752 are shown by filled circles. 



In order to compare the present models with previous 
ones (e.g., Fenner et aL 2004, Ventura & D'Antona 2005c), 
we compare the present results with the same observational 
ones (Yong et al. 2003 for NGC 6752) that have been al- 
ready compared also with these previous studies in §3. We 
also compare the present results with other more recent ob- 
servations on abundances of stars (in particular, less evolved 
stars) in GCs (e.g. Ramirez & Cohen 2002; Briley et al. 
2004a,b; Carretta et al. 2005; Cohen et al. 2005) and assess 
the viability of the external pollution scenario later in §4.1 



3 RESULTS 

3.1 The fiducial model 

3.1.1 C and N abundances 

Fig. 2 shows the time evolution of the fiducial model Ml in 
comparison with the model M2 on the [C/Fe]-[N/Fe] plane. 
There are differences in initial abundances of infalling gas 
between models Ml and M2 (e.g., [C/Fe] and [N/Fe]). Since 
these abundances are estimated for new stars formed at each 
time step in a GC, each location of the line represents abun- 
dances of stars with given ages rather than the mean abun- 
dances of the GC. As the proto-GC cloud is polluted by 
ejecta from more massive field AGB stars (mi > 5M0), 
[N/Fe] and [C/Fe] in the cloud rapidly increases and de- 
creases, respectively. Gas infall increases the total gas mass 
of the cloud (and thus the density) and consequently in- 
creases gradually star formation in the cloud. As a result of 
this, new stars that formed from the gas, which is a mix of 
infalling gas and AGB ejecta, can have higher [N/Fe] and 
lower [C/Fe]. As a larger amount of AGB ejecta is accumu- 
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Figure 3. Normalized number distributions of [C/Fe] (upper) 
and [N/Fe] (lower) in the fiducial model Ml. The [C/Fe] distri- 
butions are shown separately for N-poor stars with [N/Fe] < 1 
(blue) and for N-rich ones with [N/Fe] > 1 (red) so that a clear 
diff'erence in the [C/Fe] distributions between the two populations 
can be seen. 



lated in the cloud at later times, new stars formed later can 
have higher [N/Fe] and lower [C/Fe] than those formed ear- 
lier. Thus stars with higher [N/Fe] can show lower [C/Fe] 
in the fiducial model Ml and this anticorrelation between 
[N/Fe] and [C/Fe] can be seen in most models. 

Although the fiducial model can qualitatively explain 
the observed anticorrelation in Fig. 2 (e.g., Norris et al. 1981; 
Smith & Norris 1993), it does not pass the observed location 
of the N-rich population in NGC 6752. This quantitative in- 
consistency between the model and the observation can be 
seen in the comparative model M2 with initially small [N/Fe] 
in the infalling gas that can significantly dilute [N /Fe] during 
star formation within GCs. In the present model, [N/Fe] al- 
ways becomes significantly high (~ 1) when [C/Fe] becomes 
low irrespective of the infall rate of N-poor gas. Since only 
two data points (i.e., two stars) shown in Fig. 2 are obviously 
not enough to make any robust conclusions on the quanti- 
tative consistency of the present model with observations, 
we discuss the consistency of the model by using more data 
sets later in §4.1. 

As a result of gradual increase (decrease) of [N/Fe] 
([C/Fe]) due to AGB pollution, new stars formed at dif- 
ferent times can have different [N/Fe] and [C/Fe] and thus 
show abundance inhomogeneity in [N/Fe] and [C/Fe]. Fig. 3, 
showing the normalized number distributions of stars with 
different [N/Fe] and [C/Fe] in Ml, demonstrates that the 
[N/Fe] distribution has a clearer bimodality than the [C/Fe] 
distribution. Most N-rich stars with [N/Fe] > 1 show a 
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Figure 4. The same as Figure 2 but for Y (helium abundance 
by mass) and [N/Fe] for Ml. 
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Figure 5. The distribution of Y (hchum abundance) norniahzcd 
to the maximum values (referred to as -Fx in what follows) at the 
time T = 1 Myr (solid) and 10 Myr (dotted) for Ml. 



lower [C/Fe] (< —0.3), as expected from chemical pollution 
from massive AGB stars. Owing to the adopted density- 
dependent star formation rate, the gas infall causes two 
peaks in the star formation rate in Ml. New stars formed 
during the first peak of star formation can have lower [N/Fe] 
whereas those formed during the second peak can have 
higher [N/Fe] due to a higher degree of AGB pollution. The 
reason for the apparent absence of the bimodality in [C/Fe] 
is that [C/Fe] can not dramatically change within a short 
time scale (~ lO'^ yr) so that new stars can show a much 
narrower [C/Fc] distribution (i.e., the bimodality can not 
be evident for the adopted binning). The final mean [C/Fe] 
and [N/Fe] at T = 10^ yr are -0.09 and 0.85, respectively. 



3.1.2 He abundance 

As the gas mass of the proto-GC cloud increases owing to ac- 
cumulation of AGB material in Ml, the He abundance (Y) of 
the gas also increases in a rnonotonic fashion: the final mean 
Y at T = 10^ yr is 0.26 whereas the final Y of stars formed 
atT = 10^ yr is 0.33. As a result of this, new stars formed 
later will have higher Y than those formed earlier. Since 
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Figure 6. The same as Figure 4 but for [C/Fc]-Y (top), [N/Fc]-Y 
(second from top), [0/Fe]-Y (second from bottom), and [Na/Fe]- 
Y (bottom) for Ml. 



AGB cjocta increases [N/Fc] rapidly, new stars that form 
later and thus have higher Y can show higher [N/Fe]. Fig. 4 
clearly shows a positive yet non-linear correlation between 

Y and [N/Fe] in new stars of GCs for Ml. This correlation 
can be seen in almost all models in the present study and 
thus is simply a result of H-burning. 

Owing to the time-evolution of Y in the gas cloud, new 
stars formed at different times show different Y, which re- 
sults in a Y inhomogeneity in the GC. Fig. 5 clearly shows 
that as the proto-GC cloud is polluted by AGB ejecta to 
a larger extent, newly formed stars show a larger degree of 

Y inhomogeneity. Furthermore, the location of the peak in 
the Y distribution is also shifted to the higher values as the 
AGB pollution proceeds: the peak value changes from 0.23 
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Figure 7. The same as Figure 2 but for the [Mg/Fe]-[A1/Fe] 
relation for Ml. For comparison, observations from Yong et al. 
(2003) are shown by filled circles. The original values of abun- 
dances listed in Yong et al. (2003) are used for the observational 
plots. 
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Figure 9. The dependences of [C/Fe] (left four) and [N/Fe] 
(right four) distributions (Fx) on s. 
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Figure 8. The same as Figure 2 but for the [0/Fe]-[Na/Fe] 
relation for Ml. For comparison, observations from Yong et al. 
(2003) are shown by filled circles. 

at T = 1 Myr to 0.25 at T = 10 Myr. These results imply 
that the time scale of GC formation is an important factor 
which can determine Y distributions of GCs. 

Since abundances in the present model either increase or 
decrease in a monotonic fashion as AGB pollution proceeds, 
strong correlations or anticorrelations can be expected be- 
tween different abundances of stars in a GC. Fig. 6 shows 
that Y can anticorrelate with [C/Fe], [O/Fe], and [Na/Fe], 



though the dependences are not necessarily linear. These 
correlations can be expected for any chemical evolution 
models of GCs using stellar yields of AGB stars with no 
3rd dredge-up. 

3.1.3 Mg-Al and Na-0 anticorrelations 

Fig. 7 shows the time evolution of Ml in the [Mg/Fc]-[A1/Fe] 
plane and the observations of bright red giant stars in NGC 
6752 (Yong et al. 2003). As a natural result of ongoing chem- 
ical pollution of the proto-GC cloud by field AGB stars, the 
[Mg/Fe] and [Al/Fc] ratios of new stars within the GC de- 
crease and increase, respectively, as AGB pollution proceeds. 
Although the observed [Mg/Fe]-[A1/Fe] anticorrelation can 
be well reproduced qualitatively, the observed slope of the 
anticorrelation and the ranges of [Mg/Fe] and [Al/Fe] are not 
consistent with the predicted values: a shallower slope of the 
anticorrelation, a wider range of [Mg/Fe], a smaller range of 
[Al/Fe] axe clearly seen in the model. This inconsistency can 
be seen in the model M3 with initially higher [Al/Fe]. Al- 
though all models in the present study show [Mg/Fc]- [Al/Fe] 
anticorrelations, none of the models shows slopes as steep as 
the observed one for NGC 6752. 

Fig. 8 shows the time evolution of Ml in the [O/Fe]- 
[Na/Fe] plane as well as the observational plots of bright 
red giant stars in NGC 6752. The model predicts a positive 
correlation between [O/Fe]- [Na/Fe], which is obviously in- 
consistent with the observed [O/Fe]- [Na/Fe] anticorrelation. 
Furthermore, the predicted [Na/Fe] is too low to match the 
observed values, and this inconsistency can be seen in all 
models of the present study. This suggests that the derived 
inconsistency is due to the adopted AGB yield models in 
which no third dredge-up is assumed. Also the present study 
with [O/Fe]- [Na/Fe] correlation is in contrast with Fenner 
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Figure 10. The dependences of Y distributions (Fx) on s. 

et al. (2004) in which third dredge-up is assumed and the 

derived [0/Fo]-[Na/Fe] anticorrelation is not quantitatively 
consistent with the observed one by Yong et al. (2003) . The 
fact that both Fenner et al. (2004) and the present study 
fail to explain the observed [0/Fe]-[Na/Fe] anticorrelation 
quantitatively may well suggest that AGB pollution scenar- 
ios have a serious problem in terms of explaining the ob- 
served abundance inhomogeneity in GCs. We will discuss 
this point later in §4.1. 

3.2 Parameter dependences 

3.2.1 The dependence on s 

Fig. 9 shows the dependences of [C/Fe] and [N/Fe] number 
distributions normalized to their maximum values (denoted 
as Fx) on s, i.e., the ratio of Mm (the total mass of infalling 
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Figure 11. The same as Figure 9 but as a function of ai. 



gas) to Magb (the total mass of AGB ejecta) and thus con- 
trols the degree of AGB pollution. Fig. 9 clearly shows that 
both the shapes and the peak values of the distributions de- 
pend on s, which suggests that s is an important parameter 
for [C/Fc] and [N/Fe] distributions. Bimodal distributions 
can be secu iu [N/Fe] for models with 3 < s < 24 and the 
fraction of N-rich stars with [N/Fe] > 1 is higher in the mod- 
els with smaller s (i.e., larger mass fraction of AGB ejecta). 
As s becomes larger, both the difference between smaller 
and larger peaks in [N/Fe] and Fx in the smaller peak be- 
come larger. The derived diverse distributions of [N/Fe] can 
be compared with the observed diversity in the CN distri- 
butions of GCs (Norris 1988). 

The [C/Fe] distributions, on the other hand, do not 
depend so strongly on s and the four models show almost 
monomodal distributions. This reflects the fact that chemi- 
cal pollution by AGB stars does not change [C/Fe] as dra- 
matically as [N/Fe]. The spread in [C/Fe] is however depen- 
dent on s in the sense that the spread is larger for mod- 
els with smaller s (i.e., a larger degree of AGB pollution). 
The predicted bimodal [N/Fe] distribution and monomodal 
[C/Fe] one are discussed later in §4.1 by comparing these 
predictions with the latest observations for less-evolved stars 
(e.g. MS stars). 

Fig. 10 shows that Y distributions normalized to the 
maximum values (Ft) depend on s such that both the peak 
values and the spreads in the Y distributions are larger for 
the models with smaller s (a larger degree of AGB pol- 
lution). These results accordingly suggest that GCs with 
larger degrees of Y inhomogeneity will show higher Y. The 
results shown in Figs. 9 and 10 therefore suggest that GCs 
with higher fractions of N-rich stars with [N/Fc] > 1 can 
show higher Y and larger degrees of Y inhomogeneity. Since 
the fraction of N-rich stars in a GC can be observationally 
quantified by an "r" parameter by counting the number frac- 
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Figure 12. The same as Figure 10 but as a function of ai. 



tion of CN-strong stars vs CN-weak stars (Norris 1988), the 
above predicted correlation can be tested by investigating 
correlations between the r parameters and observable prop- 
erties dependent on Y inhomogeneity, such as color spreads 
between MS stars in the color-magnitude diagrams of GCs. 



3.2.2 The dependence on a\ 

Fig. 11 shows that [N/Fe] distributions depend strongly on 
(Ti, which means that [N/Fe] distributions can be controlled 
by the time-scale of gas infall. As seen in the s dependences 
of [N/Fc], the models with higher peak values of [N/Fe] show 
larger degrees of [N/Fe] spreads in the distributions. It is 
significant that the model with ai = 1.0 shows an almost 
monomodal [N/Fe] distribution with an apparently small 
dispersion of [N/Fe] (i.e. no clear bimodality). [C/Fe] distri- 



butions do not depend so strongly on ai as [N/Fe] so that 
models show a relatively narrow range of [C/Fej. It can be 
concluded from Figs. 9 and 11 that (i) [C/Fe] distributions 
are highly likely to be monomodal irrespective of s and cri 
and (ii) [N/Fe] distributions are likely to be bimodal and it 
depends on s and ai how clearly the distributions show the 
bimodality. 

Fig. 12 shows that the models with smaller cri (i.e., more 
rapid gas infall) show larger degrees of Y inhomogeneity. 
However a positive correlation between the peak Y and the 
degrees of Y inhomogeneity seen in models with different s 
(shown in Fig. 10) can not be seen in these models with dif- 
ferent (Ti. Given the fact that mean Y should not be different 
between models with the same A/agb (or s) yet different ai, 
r parameters can be correlated only with the degrees of Y 
inhomogeneity (i.e., not with Y). The results shown in Figs. 
10 and 12 thus suggest that Y distributions are diverse de- 
pending on the two key parameters s and cti. Observational 
implications of these results are discussed later in §4.4. 

Thus the present study shows that Y, [C/Fe], and 
[N/Fe] distributions in GCs are predicted to be diverse if 
GCs are formed from the mixture of infalling gas and AGB 
ejecta of field AGB stars in forming dwarfs at high z. The 
results shown in Figs. 10-12 can provide some physical basis 
not only for the observed diversity in CN distributions of 
GCs (e.g., Norris 1988) but also for the possible Y spreads 
in NGC 2808 (e.g., D'Antona et al. 2005) and u) Cen (e.g., 
Bedin et al. 2004) . It should be stressed here that the diverse 
distributions come from the combination of the external pol- 
lution scenario and the adopted AGB yield: neither the exter- 
nal pollution scenario assuming the third dredge-up in AGB 
stars nor classical self-pollution scenarios assuming no third 
dredge-up are able to predict diverse distributions similar to 
the observed ones. 



4 DISCUSSION 

4.1 Comparison with observations 

It is more instructive for the present study to compare our 
model predictions on abundance inhomogeneity within GCs 
with observations for GC MS stars, rather than with the gi- 
ant branch populations. Comparison with near-MS stars is 
to be preferred, because we do not have to consider abun- 
dance changes due to deep mixing in evolved stars. That 
said, there are relatively few studies that show observational 
data sets on abundances of MS populations in GCs (due to 
the difficulty in observing these faint objects by the current 
generation of telescopes). We thus use available near-MS 
data sets in order to understand how well the "pure" AGB 
pollution scenario can explain the observations for Galactic 
globular clusters, and note that the data for evolved stars 
(i.e., RGB and AGB stars) are not ideal for testing our mod- 
els. We here mainly focus on relations between [C/Fe] and 
[N/Fe] using recent observational data sets from Briley et 
al. (2004a, b), Carretta et al. (2005), Cohen et al. (2002), 
Cohen et al. (2005), and Smith et al. (2005). 

Fig. 13 shows the comparison between the predicted 
evolution of model in the [C/Fe]- [N/Fe] plane and the ob- 
served [C/Fej and [N/Fe] values for the red giants in M13, 
MIO, and NGC 7006 (Smith et al. 2005). It is clear from this 
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Figure 13. Comparison between models (Ml and M15) and 
observations (Smith et al. 2005) for MIO (red triangles), M13 
(blue squares), and NGC 7006 (green pentagons) in the [C/Fe]- 
[N/Fe] plane. The time evolution of Ml and M15 are shown by 
thick solid and dotted lines, respectively. The observational data 
are for red giants in these GCs (i.e., for more evolved stars). 



figure that (i) the observed anticorrelation between [C/Fe] 
and [N/Fe] can be qualitatively reproduced, (ii) the adopted 
"pure" AGB-poUution scenario, however, can not explain 
the observed smaller [C/Fe] ([C/Fe] < —0.5), and (iii) it can 
not explain stars with high [N/Fe] (> 1.5) and low [C/Fe] 
(< —0.5). These results clearly suggest that a combination 
of deep mixing (e.g., see the new mechanism discussed by 
Eggleton et al 2006) and primordial pollution from AGB 
stars in proto-GC clouds need to be considered for reason- 
able explanations of the observations. 

Fig. 14 shows the comparison between the predicted 
evolution of the models on the [C/Fe]- [N/Fe] plane and the 
observed [C/Fe] and [N/Fe] for less evolved stars such as 
MS populations, subgiants, and dwarfs in NGC 6397 (Car- 
retta et al. 2005), NGC 6752 (Carretta et al. 2005) , 47 Tuc 
(Briley et al. 2004b; Carretta et al. 2005), M15 (Cohen et al. 
2005), M13 (Briley et al. 2QQ4a), and M5 (Cohen et al. 2002). 
Fig. 13 and Fig. 14 clearly indicate that the models are 
much more consistent with observations for less evolved stars 
(Fig. 14) than for evolved ones (Fig. 13), though the mod- 
els still have difficulty in explaining stars with low [C/Fe] 
(< —0.5). Furthermore, the observed [C/Fe]-[N/Fe] anticor- 
relation and the apparently non-linear dependence of [N/Fe] 
on [C/Fe] are both consistent with the theoretical predic- 
tions. These derived consistencies may well provide support 
for the adopted star formation histories of proto-GC clouds 
in the external pollution scenario. 

Fig. 15 shows the comparison between the predicted 
[C/Fe] and [N/Fe] distributions and the observed ones for 
subgiants and stars at the base of RGB in M15 and for MS 
stars in 47 Tuc. It is clear from Fig. 15 that (i) both mod- 
els and observations show almost monomodal distributions 
in [C/Fe], (ii) the predicted peak (in normalized number 
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Figure 14. Comparison between models (Ml and M15) and 
observations for dwarfs and subgiants in NGC 6397 (red open 
triangles; Carretta et al. 2005), for those in NGC 6752 (red filled 
triangles; Carretta et al. 2005), for those in 47 Tuc (blue open 
squares; Carretta et al. 2005), for SGB stars in M5 (blue filled 
squares; Cohen et al. 2002), for main-sequence turn-off stars in 
M13 (green open pentagons; Briley et al. 2004a), for MS stars 
in 47 Tuc (green filled pentagons; Briley et al. 2004b), and for 
subgiants and stars at the base of the red giant branch in M15 
(magenta filled circles; Cohen et al. 2005) in the [C/Fe]- [N/Fe] 
plane. The time evolution of Ml and M15 are shown by thick 
solid and dotted lines, respectively. The typical error bars are 0.2 
for [C/Fe] and [N/Fe] (see Fig. 7 in Carretta et al. 2005). 



fraction) in the [C/Fe] distribution is significantly higher 
than those of the observations, (iii) the predicted bimodal 
distribution is more similar to the observed one in 47 Tuc 
than to the one in M15, and (iv) the locations of the two 
peaks in [N/Fe] are however not exactly the same as the 
observed ones. Thus the observed wider spread in [C/Fe] 
can not be explained so well by the present models, though 
the monomodal [C/Fe] and bimodal [N/Fe] distributions are 
quite well reproduced by the models. 

The present models also predict the mean values of 
[C/Fe] and [N/Fe] in GCs and thus can be compared with 
the ones derived from high-resolution spectroscopic studies 
of integrated spectra of GCs in nearby galaxies other than the 
Galaxy. Ongoing and future observational studies based on 
large ground-based telescopes with high-dispersion spectro- 
graphs enable us to derive [C/Fe] and [N/Fe] of GCs in other 
galaxies by model fitting (e.g., Bernstein 2006). Recent pre- 
liminary results reported by Strader et al. (2006) for a GC 
in M31 (034-096) show that [C/Fe] and [N/Fe] are roughly 
dex and 0.6 dex, respectively, with the observational errors 
of about 0.15 dex. This result is more consistent with the 
present models with smaller s (that shows higher [C/Fe] and 
smaller [N/Fe]). We suggest that a statistical study of M31 
GCs in terms of [C/Fe] and [N/Fe] can provide new insight 
into the differences in chemical pollution of proto-GC clouds 
by AGB stars between M31 GCs and the Galactic ones. 
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Figure 15. Comparison of [C/Fe] and [N/Fe] distributions be- 
tween the model Ml (red solid) and observations for subgiants 
and stars at the base of the red giant branch in M15 (blue dot- 
ted; Cohen et al. 2005) and for MS stars in 47 Tuc (green dashed; 
Briley et al. 2004b). 
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Figure 16. The predicted color spreads (A(y — /)) due to Y 
inhomogeneity for MS populations in Ml at T=l Myr (solid) and 
10 Myr (dotted). A{V — I) is calculated from the results shown 
in Fig. 5 using the formula by D'Antona et al. 2005 (see the main 
text for details). This figure implies that if GCs are formed from 
proto-GC clouds chemically polluted by AGB stars, they can show 
color spreads (up to 0.05 mag in y — /) in their MS populations. 



4.2 Advantages and disadvantages of the external 
pollution scenario 

4-2.1 Advantages 

One of the serious problems in any self-pollution scenario 
is that the number fraction of N-rich, C-depleted stars in a 
GC for a canonical IMF can not be large enough to explain 
the observed fraction (e.g., Smith & Norris 1982; D'Antona 
& Caloi 2004; Prantzsos & Charbonnel 2006). This prob- 
lem can be solved if we assume peculiar top-heavy IMFs 
(D'Antona & Caloi 2004), for which GCs can either be eas- 
ily destroyed by the Galactic tidal field during their evolu- 
tion (e.g., ChernofT & Weinberg 1990) or be disintegrated 
quickly due to loss of a large amount of mass through type- 
II supernovae explosion (Bekki & Norris 2006). Since both 
N-rich and N-normal components can be formed from proto- 
GC clouds chemically polluted by field AGB stars (that do 
not become parts of GCs) in the external pollution scenario, 
there are no serious problems on number fraction of N-rich, 
C-depleted stars in the scenario with a canonical IMF. Re- 
cently Pasquini et al. (2005) have suggested that the whole 
proto-GC cloud needs to be chemically enriched by a previ- 
ous generation of stars based on the observational data of 
the Li abundance in 9 turn-off stars in NGC 6752: this is 
consistent with the external pollution scenario. 

GCs can be formed during or after the formation of the 
field stars in dwarfs at very high z in the external scenario. 
A consequence of this is that the field stars in dwarfs on av- 
erage should be older and more metal-poor than GCs when 
GCs are formed in this scenario. Accordingly, if dwarfs with 
GCs can be destroyed by the tidal field of the proto-Galaxy 
to form the Galactic field halo stars and the halo GCs soon 
after GC formation, the field halo stars are highly likely to 
be more metal-poor than the halo GCs in this scenario. It 
is a well known observational fact that the metallicities of 
the most metal-poor halos stars are lower than those of the 
most metal-poor GCs in the Galaxy (e.g.. Freeman 1993). 
The earlier formation of field halo stars in comparison with 
GCs in the external pollution scenario can naturally explain 
this observation. 

Since proto-GC clouds are chemically enriched by old 
field stars (not by the very first stars like Pop III stars) of 
low-mass dwarf galaxies in the present models, the clouds 
can have abundances similar to field stars (which can fi- 
nally become halo stars) and be uniformly enriched by r— 
and s— process syntheses. Therefore, GCs formed from these 
clouds can have abundances similar to those of field stars 
except for light elements and accordingly explain recent ob- 
servations that have reported little differences in [Fe/H]- 
dependences of [Sr/Fe], [Y/Fe], [Ba/Fe], and [Eu/Fe] in main 
sequence turnoff stars and subgiants at the base of the Red 
Giant Branch (RGB) in some GCs (James et al. 2004). Stars 
in the models with no third dredge-up keep initial abun- 
dances of s— process elements so that the models can be 
consistent with the observed possible evidences for homo- 
geneous abundances of s-process elements in James et al. 
(2004). The external pollution from AGB stars with third 
dredge-up would be inconsistent with the above observations 
by James et al. (2004) : future observations on the degree of 
homogeneity in the abundances of s— process elements will 
give strong constrains on any AGB-poUution scenario. 

Recent abundance studies of the Galactic open clus- 
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ter have not yet revealed possible abundance spreads be- 
tween stars of these clusters (e.g., Friel et al. 2003; Dc Silva 
et al. 2006). The external pollution scenario suggests that 
these clusters are highly unlikely to show abundance inho- 
mogeneitios, because they arc formed from the Galactic gas 
clouds in which AGB pollution can not proceed owing to 
very shallow gravitational potentials of the clouds. Hill et al. 
(2000) reported that LMC clusters do not show Al-O anti- 
correlations, which implies that AGB pollution did not pro- 
ceed in these systems. The external pollution scenario thus 
suggests that these LMC clusters, which shows disky kine- 
matics (e.g.. Freeman et al. 1983), were formed not within 
dwarfs that were the LMC's building blocks but within the 
proto-LMC disk via a dissipativc process at the formation 
epoch of the disk. The observed differences in some elemen- 
tal abundances between the LMC clusters and the Galactic 
GCs (Johnson et al. 2006) may well be understood in terms 
of clusters originating within the forming disk of the LMC. 



4.2.2 Disadvantages 

As pointed out by Fenner et al. (2004), chemical evolu- 
tion models of GCs based on AGB-poUution have some 
difficulties in explaining quantitatively O-Na and Mg-Al 
anticorrelations and correlations between ■^^Mg/^'*Mg (and 
^'^Mg/^'iMg) and light elements (e.g., [0/Fe] and [Na/Fe]). 
Although some improvement over previous models can be 
seen in the present models in which AGB stars experience 
HBB but no third drcdge-up, the observed O-Na anticorrcla- 
tion can not be explained by them. This incapability comes 
from the adopted stellar yield models of AGB stars rather 
than from improper models for early star formation histories 
within GCs. Given the fact that there are some uncertain- 
ties in model calculations of yields in AGB stars (Ventura 
& D'Antona 2005c; Campbell et al 2006), the above appar- 
ent failure does not necessarily mean that there is a serious 
problem with the adopted AGB-poUution scenario. Indeed, 
we remind the reader that no consistent solution currently 
exists. All proposals rely on tweaking of some inputs, be it 
overshoot, reaction rates or something else, or indeed a com- 
bination. One solution for the above problem would be to 
adopt AGB models with little third dredge-up (i.e., not the 
models with no dredge-up) in order to avoid unreasonably 
small values of [Na/Fe] in the models. This contrived way 
of modeling, however, should not be regarded as a neces- 
sary solution for this problem. It may say something about 
the characteristics of low-metallicity AGB stars or perhaps 
higher mass super-AGB stars - maybe they do not suffer as 
much third dredge-up as expected by our models. 



4.3 The diversity in CN distributions 

Previous observational studios showed that there exists vari- 
ations in the strength of CN bands between the Galac- 
tic GCs (e.g., Norris & Smith 1981; Suntzeff 1981). Norris 
(1988) investigated the value of the "r" parameter, which is 
defined as the ratio of the number of CN-strong to CN- 
weak stars in a given system, for each of 12 GCs with 
— 1.9 < [Fe/H] < —1.2 and thereby found that the values 
of the r parameter range from 0.22 (i.e., smaller fraction of 
CN-strong stars) to 3.29 (i.e., larger fraction). Norris et al. 



(1984) also suggested that the observed anticorrclation be- 
tween the behavior of the CN bands and that of the CH fea/- 
tures in GCs (e.g., Norris & Cottrell 1979; Norris & Freeman 
1982; Norris et al. 1984; Bell et al. 1983) is difficult to ex- 
plain in terms of the primordial pollution (or self-pollution) 
scenarios. 

We have demonstrated that (1) the distributions of N 
and C are diverse depending on the two key parameters s 
and (71, (2) N distributions show clear bimodality in some 
models, and (3) C and N abundances aiiticorrclato with each 
other in the external pollution scenario. These results pro- 
vide some physical basis for the origin of the observed di- 
versity and for the observed C-N anticorrclation, if we un- 
derstand what determines the values of the above two key 
parameters in proto- galaxies. Owing to the lack of extensive 
numerical simulations on the dynamical fate of AGB ejecta 
in the central region of proto-galaxies, it is unclear what 
physical processes in proto-galaxies (e.g., the effectiveness 
of stellar feedback effects) can control s and oi. We plan to 
derive s and cri and their dependences on physical properties 
of proto-galaxies (e.g., masses and sizes) in our future high- 
resolution simulations and thereby discuss the origin of the 
observed diversity of CN distributions in GC in the context 
of galaxy formation. 



4.4 Possible helium abundance inhomogeneity in 
GCs 

Recent photometric observations of stars in ui Cen have dis- 
covered a double main sequence (DMS) in the color magni- 
tude diagrams (CMDs) of its stellar populations (Anderson 
1997; Bedin et al. 2004). One of the most promising inter- 
pretation is that stars on the bluer main sequence (bMS) of 
the DMS represents a very helium-rich {Y > 0.3) population 
(Norris 2004; Piotto ct al. 2005). Furthermore, D'Antona et 
al. (2005) have reported helium abundance variation among 
main sequence stars of NGC 2808. Green (1980) first showed 
that there is a possible correlation between the helium abun- 
dance of GCs and the locations of GCs with respect to the 
Galactic center: imicr GCs arc more likely to show higher he- 
lium abundance. Recently Yong et al. (2005) have reported a 
possible trend between [Fe/H] and [Al/Fe] for giants stars in 
NGC 6752 and thus suggested that this trend can result from 
a possible correlation between the H-He ratio and [Al/Fe], 
which is a natural result of chemical pollution by H-burning 
ejecta. Although the total number of observations reporting 
the possible abundance spread in helium is currently very 
small, these observations contain valuable information not 
only on horizontal branch morphologies but also on helium 
pollution processes in proto-GCs (e.g., D'Antona et al. 2005, 
2006). 

Important predictions of the present external pollution 

scenario on helium abundance properties of GCs are (1) 
spread in helium abundance between stars, (2) larger helium 
abundances in N-rich, C-dcplctcd stars, and (3) smaller he- 
lium abundance in GCs formed from gas less polluted by 
field AGB stars. Furthermore, we can provide some predic- 
tion on color distributions of stars for the main-sequence 
(MS) populations of GCs by using the model developed by 
D'Antona ct al. (2005). Differences in Y between stars with 
normal Y (=0.24) and those with higher/lower Y can be con- 
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verted into color difference (A(V — /)) by using the following 
formula (D'Antona et al. 2005): 



A{V-I) 
AY 



-0.438. 



(9) 



Fig. 16 clearly demonstrates that the color distributions of 
stars in the fiducial model Ml shows significant dispersions 
in F — / of the MS. This dispersion can be observed as 
spreads in the CM relations of MS populations of GCs and 
thus as evidence for GC formation from AGB cjccta. 

Several authors have pointed out that He abundance 
inhomogeneity of stars in GCs can be imprinted on pho- 
tometric properties of stellar populations on the horizontal 
branch (HB) of GCs (e.g., Norris 2004; D'Antona et al. 2006; 
Lcc ct al. 2005). Those authors suggested that a GC with 
a He-rich population can (e.g., NGC 2808) show a charac- 
teristic HB color and luminosity distribution, a bimodal HB 
distribution, and faint, very hot HB stars. Since a GC with 
a larger fraction of N-rich stars can have a higher degree of 
He abundance inhomogeneity in the external pollution sce- 
nario, the HB properties of GCs can be correlated with the 
fraction of N-rich stars (i.e., r parameter). If the origin of 
the second parameter problem is related to the difference 
in helium abundance between GCs (see Catlan 2005 for a 
recent review), the external pollution scenario suggests that 
the origin of the problem can be closely associated with dif- 
ferences in AGB pollution processes (controlled mainly by 
the two parameters, s and ai) between different proto-GC 
clouds in the central regions of forming dwarfs at high z. 



4.5 Origin of the observed correlation between 

the number fraction of CN-strong populations 
and ellipticities of GCs 

Norris (1987) investigated possible correlations between the 
r parameter and other physical properties of GCs (e.g., 
metallicities and luminosities) and found that there exists 
a correlation between ellipticities (e) of GC shapes and the 
values of the r parameters (See Smith (2002) for recent con- 
firmation of this correlation). This result means that if the 
flattened shapes arc due to rotational kinematics of GCs 
(rather than anisotropy of velocity dispersion), GCs with 
a larger amount of angular momentum are more likely to 
show larger r. Norris (1987) suggested that if global rota- 
tion of a GC positively correlate with rotation of individual 
stars within the GC and if evolutionary mixing responsible 
for the GN abundance inhomogeneity is driven by stellar ro- 
tation, the observed correlation can be readily explained. 
Norris (1987) however pointed out that this explanation has 
a difficulty in explaining the existence of CN-enhanced stars 
near the main-sequence turnoff of 47 Tuc (Bell et al. 1983) 
unless the effects of rotation are efficient for main-sequence 
stars. Furthermore it is unclear why stars within GCs with 
a larger amount of intrinsic angular momentum would have 
a larger amount of stellar rotation in this scenario by Norris 
(1987). 

In the present external pollution scenario, the values of 
the r parameters depend strongly on the relative mass frac- 
tion of AGB ejecta with respect to the total gas mass (i.e., 
AGB ejecta -|- fresh infall gas) in proto-GCs: The higher the 
fraction is, the larger the r becomes. One of possible expla- 
nations for the above e-r correlation in the context of the ex- 



ternal pollution is therefore that whenever a larger amount 
of AGB ejecta can be transferred into proto-GC regions (i.e., 
central regions of their host dwarfs), the gas is highly likely 
to have a larger amount of angular momentum, from which 
GCs with a larger amount of angular momentum and thus 
with larger e can be formed. We suggest that the above sit- 
uation of the larger amount of AGB ejecta with higher (spe- 
cific) angular momentum is likely if field stars in the central 
regions of dwarfs (thus, AGB polluters for proto-GCs) have 
rotational kinematics: AGB ejecta from such stars rotating 
dwarfs can also have orbital angular momentum (with re- 
spect to the centers of dwarfs) so that GCs formed from 
the ejecta in the centers can have rotational kinematics and 
thus large e after conversion of orbital angular momentum 
of gas into intrinsic spin angular momentum of GCs during 
GC formation. 

It is thus our future study to investigate (i) whether 
orbital angular momentum of AGB ejecta from field stars 
rotating dwarfs can be really converted into intrinsic spin an- 
gular momentum of GCs and (ii) whether GCs have a larger 
amount of intrinsic spin angular momentum thus larger e, 
if they are formed from gas with a larger fraction of AGB 
ejecta. High-resolution numerical simulations on stellar and 
gas dynamics in the central 0.1 — lOOpc of dwarfs at high 
z would be essential for this investigation. We plan to per- 
form such chemodynamical simulations with AGB feedback 
effects (e.g., thermal heating and return of AGB ejecta to 
ISM) to understand the origin of the observed e-r correla- 
tion in the Galactic GCs. Nuclear dynamics in proto-galaxies 
would be dependent on the depth of their gravitational po- 
tential wells thus on their masses and sizes. Our future sim- 
ulations will therefore enable us to understand whether very 
flattened GCs (e.g., uj Ccu) can be formed in less or more 
massive dwarfs at high z. 

4.6 Origin of possible higher nitrogen abundances 
for some bright GCs in M31 

Recently, Li & Burstein (2003) have found that the NH ab- 
sorption line is far stronger in (three) metal-rich GCs of 
M31 than it is in the Galactic GCs at a given value of CH or 
[Fe/H]. Beasley et al. (2004) also have revealed that the near- 
UV cyanogen features of M31 GCs are strongly enhanced 
with respect to the Galactic GCs for -1.5 < [Fe/H] < -0.3. 
These observations imply that [N/Fe] can be higher in some 
M31 GCs than in the Galactic ones for a same metallic- 
ity range. It has been also revealed that M31 GCs hav- 
ing very strong NH absorptions are more luminous with 
-11 < My < -8.5 (e.g., Burstein et al. 2004). This implies 
that the possibly higher [N/Fe] is due to a larger rmmber 
of bright GCs in M31, which can have the higher values for 
some physical reasons. It remains, however, totally unclear 
why these bright M31 GCs have high [N/Fe]. 

The present models have shown that mean values of 
[N/Fe] of GCs, which should be compared with the [N/Fe] 
estimated from the integrated spectra of M31 GCs, can be 
as high as 0.9. This result, combined with recent theoretical 
results and observational ones on nucleus formation of galax- 
ies (e.g., Bckki 2006; Bckki et al. 2006; Cote et al. 2006), can 
provide the following answer for the above question. Based 
on analytical calculations, Bekki (2006) pointed out that 
AGB ejecta can be more effectively trapped in nuclear re- 
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gion of more massive dwarfs embedded in dark matter halos 
so that AGB pollution can proceeds to a greater extent in 
these dwarfs. Inward radial gas-transfer into nuclear regions, 
which is indispensable for the formation of compact stellar 
systems, is demonstrated to be more efficient in more mas- 
sive dwarfs (Bekki et al. 2006). Photometric studies on nu- 
clear properties of dwarfs based on the ACS Virgo Cluster 
Survey have recently revealed that stellar nuclei are more 
massive in more massive dwarfs (Cote et al. 2006). There- 
fore we can claim that if some bright GCs in M31 originate 
from nuclei of brighter nucleated dwarfs, the observed higher 
[N/Fe] can be naturally explained. 

Thus the origin of the observed possible differences in 
[N/Fe] between M31 GCs and the Galactic ones could re- 
flect the fact that the proto-M31 contained a larger number 
of brighter dwarfs where AGB pollution could proceed more 
effectively in their central regions to form GCs with higher 
[N/Fe]. These more massive dwarfs were destroyed to form 
the MSl's stellar halo during the hierarchical formation of 
M31 through merging of these dwarfs. Because of the mass- 
metallicity relation of field stellar populations of dwarfs (e.g. 
Durrell et al. 1996), the developed stellar halo can show a 
higher metallicity and is consistent with observations (e.g., 
Reitzel et al. 1998). If [N/Fe] of GCs can be controlled by 
nuclear dynamics of their host dwarfs at very high z, as sug- 
gested by the present study, more luminous galaxies would 
have both a larger number of GCs with high [N/Fe] and more 
mctal-rich (thus redder) stellar halos. We thus suggest that 
future high-resolution spectroscopic studies for GC abun- 
dances in galaxies beyond the Local Group, combined with 
wide-field imaging of stellar halos in these galaxies, can give 
some constraints on the formation histories of GCs within 
dwarfs at very high z. 



5 CONCLUSIONS 

We have investigated the abundance inhomogeneity among 
the light elements (e.g., C, N, O, Na, and Al) in stars in 
the Galactic GCs by adopting the external pollution sce- 
nario and using latest stellar yield models of metal-poor 
AGB stars with and without third dredge-up. In the exter- 
nal pollution scenario, GCs within a forming low-mass dwarf 
embedded in a dark matter halo can be formed from the 
mixture of (i) gas ejected from the field AGB stars formed 
earlier in the dwarf and (ii) the interstellar gas infalling to 
the central region of the dwarf. We extensively investigated 
our models with no third dredge-up over wide parameter 
ranges. We summarise our principal results of the models as 
follows. 

(1) The ratio of the total mass of infalling gas to that 
of AGB ejecta during GC formation within a forming dwarf 
galaxy (described by the s parameter) and the time scale of 
gas infall (ai) are the most important parameters that can 
determine abundance properties of GCs. 

(2) Both [N/Fe] and [C/Fe] can be diverse among stars 
within a GC owing to chemical pollution from field AGB 
stars. [N/Fe] distributions in some GCs can clearly show 
bimodality wlicrcs [C/Fe] can be monomodal in most mod- 
els. [N/Fe] distributions depend on s, such that models with 
smaller s (i.e., larger mass fraction of AGB ejecta used for 
GC formation) show the [N/Fe] bimodality more clearly. 



(3) N-rich, C-poor stars in GCs also have higher He 
abundances owing to pollution from massive AGB stars with 
He-rich ejecta. The number fraction of He-rich stars (Y > 
0.30) is higher for the models with smaller s and shorter ai 
for 3 < s < 24 and 10''' < cri < lO'^ yr. He abundances of 
stars can correlate with [N/Fe] and [Al/Fe] and anticorrelate 
with [C/Fe], [O/Fe], and [Na/Fe] within GCs (and a simple 
consequence of H-burning). 

(4) Although our models can much better explain the 
observed C-N and Mg-Al anticorrelations than previous the- 
oretical models, the observed O-Na anticorrelation can not 
by simply explained by any models with different param- 
eters in this scenario. The inability of the present models 
to match the observations results from the adopted AGB 
yields and thus suggests that models using different AGB 
yields need to be explored to assess the viability of AGB 
pollution scenarios (For more details, see Appendix B). 

(5) Extensive comparison between the models and the 
observations demonstrates that the predicted [C/Fc]-[N/Fe] 
anticorrelation can be more quantitatively consistent with 
observations for less evolved stars (MS populations, sub- 
giants, and dwarfs) than for evolved ones (e.g., red giants). 
The comparison also demonstrates that the observed [C/Fc] 
and [N/Fe] distributions for less evolved stars are at least 
qualitatively consistent with the model predictions. 

(6) One important prediction of our models is that GCs 
with higher degrees of abundance inhomogeneity should 
show a larger color spread (e.g., A(V — /)) in their MS pop- 
ulations as the result of bluer stars formed from He-rich gas. 
Accordingly, it would be worthwhile for future observational 
studies to investigate correlations between the r parameters 
(i.e., the number ratio of CN-strong to CN-weak stars) and 
A{V - I) in GCs. 

(7) The values of the two key parameters (s and ai) are 
suggested to be controlled by stellar and gas dynamics in the 
central regions of host galaxies of GCs. Central stellar and 
gas dynamics (e.g., gas fueling rates to the central 10 pc) 
have been suggested to be controlled by global properties 
(e.g., total masses and densities) of galaxies. Accordingly, 
abundance inhomogeneities of GCs may have fossil infor- 
mation on (i) formation and evolution histories of central 
regions of defunct dwarfs and (ii) physical properties of the 
dwarfs. 

Thus we have shown that the CN-bimodality and the 
number fraction of C-depleted, N-rich populations of a GC 
can give strong constants on the short-term (an order of lO'^ 
yr) star formation history of the proto-GC cloud of the GC, 
if the cloud is chemically polluted by AGB stars. We ac- 
cordingly suggest that any theoretical models for the origin 
of GC abundances need to explain not only the O-Na and 
Mg-Al anticorrelations but also the bimodal CN abundance 
distributions in GCs. 
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APPENDIX A: A POSSIBLE MODEL FOR THE 
BMS IN w CEN 

Bcdin ct al. (2004) suggested that stars on the bluer main 
sequence (bMS) of the double main sequence (DMS) can be 
a very helium-rich {Y > 0.3) population. A number of re- 
cent investigations have suggested that a higher degree of 
He enhancement in the bMS is a viable explanation for the 
origin of the DMS (e.g.. Loo ct al. 2005; Norris 2004; Pi- 
otto et al. 2005) and recent analytical works by Bekki & 
Norris (2006) have suggested that He enhancement due to 
self-pollution (or more correctly, external pollution) of AGB 
stars can explain the origin of the bMS. It however remains 
unclear whether such AGB pollution scenarios can quan- 
titatively explain the observed abundance properties (e.g., 
[C/Fe], [N/Fe], and Y) in a self-consistent manner. 

Fig. Al shows Y, [C/Fo], and [N/Fo] distributions of 
the model M14 in which very small s (i.e., a higher degree 
of AGB pollution), smaller ai (i.e., more rapid gas infall), 
and smaller Mg(0) (i.e., more rapid chemical enrichment) 
are adopted so that the model shows a large fraction of He- 
rich stars. Owing to the modeled very larger degree of AGB 
pollution, stars in this model shows a large fraction of Hc- 
rich (> 0.3) and N-rich (> 1) stars. Piotto et al. (2005) have 
reported that (i) the He-rich bMS shows a rather large nitro- 
gen abundance ([N/M] ^ 1 — 1.5) and (ii) it is not very rich 
in C ([C/M]=0). The result in Fig. Al implies that these ob- 
servations can be due to AGB pollution, though the C abun- 
dance of the bMS is not so consistent with the predicted one. 
An important result in Fig. Al is that some fraction (~ 0.1) 
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Figure Al. A possible model of Y (top), [C/Fe] (middle), and 
[N/Fe] (bottom) distributions for stars on the blue main sequence 
(bMS) in u! Cen which is suggested to have a significant fraction 
of He-rich stars in the bMS population. 



of stars show relatively normal Y (0.24 — 0.25). This sug- 
gest that if the bMS is formed from gas polluted by AGB 
stars, stars with relatively normal Y arc also formed, though 
the mass fraction of the He-normal population is small. This 
may well provide a clue to the origin of the observed nor- 
mal Y (=0.246) for the metal-intermediate RRL stars (that 
would have been formed with the bMS) in to Cen (Sollima 
et al. 2006). 

Very luminous GCs such as ui Cen in the Galaxy and 
Gl in M31 have long been suggested to originate from nuclei 
of ancient nucleated systems, where continuous gas infall 
and chemical enrichment were highly likely to occur (e.g., 
Hilkor & Richtler 2000; Bekki & Freeman 2003; Bekki & 
Chiba 2004). As demonstrated above, the total mass of AGB 
ejecta used for the He-rich population in u) Cen needs to 
be ~ 6 times larger than that of infalling gas not enriched 
by AGB stars. It remains unclear whether this can really 
happen in the nuclei of dwarf galaxies that are considered 
to be progenitors of very luminous GCs like uj Cen and Gl. 
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APPENDIX B: COMMENTS ON AGB MODELS 
WITH AND WITHOUT THIRD DREDGE-UP 

Current standard models of intermediate-mass AGB stars 
with third dredge-up (3DUP) do not reproduce the ob- 
served C-N, O-Na or Mg-Al inverse-correlations observed 
in most GCs. We note however that the models of Ventura 
& DAntona (2005b) do come close, due to their use and 
development of an alternative theory of convection, which 
alters the evolution of their models significantly (see e.g., 
Ventura & DAntona 2005a). In the standard models the 
combination of hot bottom burning (HBB) and 3DUP dur- 
ing AGB evolution does predict enhancements of nitrogen, 
sodium and aluminium, and a reduction in oxygen - all of 
which arc needed to explain the observations - however they 
also predict an increase of carbon and magnesium, giving 
positive C-N and Mg-Al correlations which is in direct vi- 
olation of the observations. In addition to this the models 
do not quantitatively account for the large degree of O de- 
pletion and tend to over-produce Na, thus falling short of 
the mark with the O-Na inverse-correlation as well (e.g., 
Denissenkov et al. 1997; Denissenkov & Herwig 2003; Ven- 
tura ct al. 2004). Furthermore the standard AGB models 
produce ejecta with an overall increase in CNO nuclei of ^ 
1 to 3 dex (mainly due to the dredge-up of carbon), whilst 
observations show that the sum of CNO nuclei is virtually 
constant, varying by a factor of only two or so at most (e.g., 
Pilachowski 1988; Smith et al. 1996; Ivans et al. 1999). These 
discrepancies are based on confronting single stellar models 
with observations. Taking the evolution of an entire cluster 
of stars into account - in particular convolving with an IMF 
- presents further problems as the stars with the prerequi- 
site hot hydrogen burning arc few in number when using a 
standard IMF (e.g., Fenner et al. 2004), whilst these stars 
need to contribute an amount of polluting material that is a 
large fraction of the GC mass. Despite all these discrepan- 
cies it should be stressed that the AGB pollution scenario 
should not be dismissed outright just yet, as the models 
are well known to have many serious uncertainties (e.g., the 
treatment of convection, overshoot, mass loss and the un- 
certainties in reaction rates) . For example the non-standard 
models of Ventura & DAntona (2005b), which use a differ- 
ent convection model, can reproduce most of the abundance 
anomalies, although some overshooting and the adoption of 
a high mass-loss rate is required. In addition to the uncer- 
tainties in the input physics of all stellar models, numerical 
details can also significantly affect the results. For example, 
time-stepping during diffusive mixing can cause changes in 
the stellar structure, hence altering the temperature at the 
bottom of the convective envelope, and thus the resulting 
nucleosynthesis and yields (Siess, private communication). 
Attempts have been made to reconcile the (standard) mod- 
els with the observations but have universally ended with a 
need to 'tweak' the relative effects of HBB and 3DUP to an 
unpalatable degree (e.g., Denissenkov & Herwig 2003), also 
suggested by Fenner et al. (2004). The models using a differ- 
ent convection formalism (Ventura & DAntona 2005b) show 
a better agreement with observations but, as mentioned ear- 
lier, also need some overshooting and the adoption of a high 
mass-loss rate. 

With the results of all these stellar modelling studies 
in mind we have chosen to investigate the limiting case by 



utilising yields from AGB models that do not experience 
any 3DUP, which we have calculated ourselves. Although 
'turning off' 3DUP in our models is ad-hoc, it was surpris- 
ingly easy to achieve. By adding an extra physical condition 
for the determination of convective boundaries (the molecu- 
lar weight gradient) via implementing the Ledoux criterion 
rather than the Schwarzschild criterion we found that no 
dredge-up at all occurred in our models. We note that we 
have not extended the convective boundaries past the formal 
ones given by the Ledoux criterion. Wo leave the detailed de- 
scription of the models for a separate paper, but now discuss 
the nucleosynthetic results relevant for the current problem. 

Decoupling 3DUP from HBB exposes the dependence 
of the yields on the interplay between the two processes and 
provides a possible solution to the mismatch between the- 
ory and observation. At the very least, by separating the 
effects it is a useful investigative tool. A salient result from 
our new stellar models, evident in Table 1, is that oxygen is 
now heavily depleted in the more massive stars. Also of im- 
portance is the fact that there is now an opposite trend for 
Mg (in the same stars) - it is depleted instead of produced. 
Furthermore, the C is depleted whilst the N yield remains 
high. All of these features derive from the lack of 3DUP, 
as O is now not periodically replenished, C can not increase 
and fresh fuel for the MgAl cycle is not available. The lack of 
3DUP also implies that the sum of CNO nuclei will remain 
constant, as fresh CNO is not added to the envelope and the 
CNO burning cycles (occurring at the base of the convective 
envelope) conserve the number of CNO nuclei. All these fea- 
tures indicate that a better fit to the observations is expected 
(IMF withstanding). Also, as stated above, these models 
may simulate the effects of a previously ignored population 
of stars, the Super- AGB stars. The main discrepancy we now 
have - in terms of comparing at a single mass - is with Na. 
At the high masses (ie. temperatures) required to deplete O 
and Mg via the MgAl and ON cycles, Na is not produced 
enough, or rather it is destroyed too much by proton cap- 
tures, resulting in low yields. Moving to lower masses (tem- 
peratures) improves the Na situation but prevents O and Mg 
from being depleted. Interestingly this is the identical prob- 
lem that Ventura & D'Antona (2005b) have with their best- 
fit (for the Na-O inverse-correlation) model. Denissenkov & 
Herwig (2003) also report this problem. We note that the 
models of Ventura & D'Antona (2005b) are quite indepen- 
dent as they use a different treatment for convection. In a re- 
cent and very timely study Ventura & D'Antona (2006) have 
investigated this problem in detail, finding that their stel- 
lar model requires either the addition of some overshooting 
or a reduction of the ^^Na(p, a)^°Ne reaction rate in order 
to produce the required Na yield whilst maintaining the O 
destruction. The Ventura and D'Antona (2005b, 2006) mod- 
els, like our more ad-hoc models, otherwise match the abun- 
dance anomalies fairly well, with the strong exception of 
the Mg isotopes, which seem resistant to any approach tried 
thus far. Note that we do not suggest that real AGB models 
of low metallicity do not suffer from the 3DUP. Rather we 
have used this as an investigative tool to simulate reduced 
3DUP, or SAGB stars with little or no 3DUP, or some other 
source of (almost) pure H-burning. Clearly a thorough un- 
derstanding, and reliable quantitative prediction of 3DUP is 
still required, for stars of all masses and metallicities. 



